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During the last 70 years, the quantum theory of angular momentum has been successfully applied 
to describing the properties of nuclei, atoms, and molecules, their interactions with each other as well 
as with external fields. Due to the properties of quantum rotations, the angular momentum algebra 
can be of tremendous complexity even for a few interacting particles, such as valence electrons 
of an atom, not to mention larger many-particle systems. In this work, we study an example of 
the latter: a rotating quantum impurity coupled to a many-body bosonic bath. In the regime 
of strong impurity-bath couplings the problem involves addition of an infinite number of angular 
momenta which renders it intractable using currently available techniques. Here, we introduce a 
novel canonical transformation which allows to eliminate the complex angular momentum algebra 
from such a class of many-body problems. In addition, the transformation exposes the problem’s 
constants of motion, and renders it solvable exactly in the limit of a slowly-rotating impurity. 

We exemplify the technique by showing that there exists a critical rotational speed at which the 
impurity suddenly acquires one quantum of angular momentum from the many-particle bath. Such 
an instability is accompanied by the deformation of the phonon density in the frame rotating along 
with the impurity. 


I. INTRODUCTION 

An important part of modern condensed matter 
physics deals with so-called ‘impurity problems’, aim¬ 
ing to understand the behavior of individual quantum 
particles coupled to a complex many-body environment. 
The interest in quantum impurities goes back to the clas¬ 
sic works of Landau, Pekar, Frohlich, and Feynman, who 
showed that propagation of electrons in crystals is largely 
affected by the quantum field of lattice excitations and 
can be rationalized by introducing the quasiparticle con¬ 
cept of the polaron mm- In turn, the properties of a 
quantum many-body system can be drastically modified 
by the presence of impurities. The most known exam¬ 
ples are the Kondo effect [5] - suppression of electron 
transport due to magnetic impurities in metals - and 
the Anderson orthogonality catastrophe which leads to 
the edge singularities in the X-ray absorption spectra of 
metals [6]. 

In many instances, the impurities - even those possess¬ 
ing an internal structure - can be accurately described as 
point-like particles. The latter is justified by the separa¬ 
tion of the energy scales inherent to the impurity and the 
surrounding bath. A well-known example is that of Bose- 
and Fermi-polarons realized in cold atomic gases by a 
number of groups There, the spherically symmet¬ 

ric ground state of an alkali atom lies hundreds of THz 
lower than any of its electronically excited states. Given 
ultracold collision energies, such an energy gap renders 
all the processes happening inside of an atom irrelevant. 

More complex systems, such as molecules, are ex¬ 
tended objects and therefore possess a number of fun- 
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damentally different types of internal motion. The lat¬ 
ter stem from the relative motion of the nuclei, such as 
rotation and vibration, which couple to each other as 
well as to the electronic spin and orbital degrees of free¬ 
dom [TTHSO]- This results in a rich low-energy dynam¬ 
ics which is highly susceptible to external perturbations. 
Moreover, in many experimental realizations molecular 
rotation is coupled to a phononic bath pertaining to the 
surrounding medium such as superfluid helium m , rare- 
gas matrix [22], or a Coulomb crystal formed in an ion 
trap [23], which needs to be properly accounted for by a 
microscopic theory. 

The concept of orbital angular momentum, however, 
goes far beyond physically rotating systems and is being 
used to describe e.g. the excited-state electrons in solids, 
whose motion is perturbed by lattice vibrations m , or 
Rydberg atoms immersed into a Bose-Einstein conden¬ 
sate P5l Despite the ubiquitous use of the angu¬ 

lar momentum concept in various branches of physics, a 
versatile theory describing the redistribution of orbital 
angular momentum in quantum many-body systems has 
not yet been developed. 

Recently, we have undertaken the first step towards 
such a theory by deriving a generic Hamiltonian which 
describes the coupling of an 5'0 (3)-symmetric impurity 
- a quantum rotor - with a bath of harmonic oscilla¬ 
tors m- We have shown that the problem can be ap¬ 
proached most naturally by introducing the quasiparticle 
concept of the 'angulon’ - a quantum rotor dressed by 
a quantum field. The angulon is an eigenstate of the 
total angular momentum of the system, which remains 
a conserved quantity in the presence of the impurity- 
bath interactions. It was found that even single-phonon 
excitations of the bath alone are capable of drastically 
modifying the rotational spectrum of the impurity, which 
manifests itself in the emerging Many-Body-Induced Fine 



Structure m 

Here we demonstrate that rotation of an anisotropic 
impurity can, in turn, substantially alter the collective 
state of a many-particle system. The effects are most 
significant in the regime of strong correlations, which 
however requires adding an infinite number of angu¬ 
lar momentum vectors pertaining to possible many-body 
states. The resulting angular momentum algebra involves 
Wigner Snj-symbols [28] of an arbitrarily high order and 
is therefore intractable using standard techniques. In or¬ 
der to overcome this problem, here we introduce a canon¬ 
ical transformation, which, to our knowledge, has never 
appeared in the literature before. The transformation 
renders the Hamiltonian independent of the impurity co¬ 
ordinates, thereby eliminating the complex angular mo¬ 
mentum algebra from the many-body problem. Further¬ 
more, the transformation singles out the conserved quan¬ 
tities of the many-body problem and renders it solvable 
exactly in the limit of a slowly rotating impurity. 

The transformation makes it apparent that there ex¬ 
ists a critical rotational speed which leads to an instabil¬ 
ity, accompanied by a discontinuity in the many-particle 
spectrum. Unlike in the vortex instability, originating 
from rotation of a condensate around a given axis [29|, 
the instability we uncover here corresponds to the fi¬ 
nite transfer of three-dimensional angular momentum be¬ 
tween the impurity and the bath. It exists solely due to 
the discrete energy spectrum inherent to quantum rota¬ 
tion. We demonstrate that the emerging instability is 
ushered by a macroscopic deformation of the surround¬ 
ing bath, i.e. the phonon density modulation in the frame 
co-rotating with the impurity. 


II. THE CANONICAL TRANSFORMATION 


The second term of Eq. § represents the kinetic en¬ 
ergy of the bosonic bath, where the corresponding cre¬ 
ation and annihilation operators, and 6k, are expressed 
in the spherical basis, bl.^^ and Here k = |k|, while 

A and /i define, respectively, the boson angular momen¬ 
tum and its projection onto the laboratory 2 ;-axis, see 
Appendix [A| for details. 

The last term of Eq. 0 describes the interaction 
between the impurity and the bath. The angular- 
momentum-dependent coupling strength, Ux{k), depends 
on the microscopic details of the two-body interaction 
between the impurity and the bosons. Eor example, in 
Ref. [27] we showed that for a linear rotor immersed into 
a Bose gas, the couplings are given by 


Ux{k)=ux 


Sk'^ekP 

Wfe(2A + 1) 


-, 1/2 


J drr'^fx{r)jx{kr). (2) 


This assumes that in the impurity frame, the interaction 
between the rotor and a bosonic atom is expanded as 


Mmp-bos(r') = Y,uxfx{r')Yxo{e',^'), (3) 

A 

with ux and /a(^ 0 giving the strength and shape of the 
potential in the corresponding angular momentum chan¬ 
nel. The prefactor of Eq. 0 depends on the bath den¬ 
sity, p, the kinetic energy of the bare atoms, e/., and 
the dispersion relation of the bosonic quasiparticles, ujk- 
Since the angulon Hamiltonian 0 describes the interac¬ 
tions between a quantum rotor and a bosonic bath of, 
in principle, any kind, we will approach it from an en¬ 
tirely general perspective, exemplifying the couplings by 
the ones of Eq. 0 . 


We start from the general Hamiltonian of the angulon 
problem, as defined in Ref. m 


H = B3'^ + ^u:J)l^J)kx^, 

kXfi 

+ ^ Ux{k) + Yx^ie, , (1) 

kXfi 


where Yx^{0, (j)) are the spherical harmonics [28] depend¬ 
ing on the molecular angle operators 0 and 0, = J dk^ 

and h= 1. 

The first term of Eq. corresponds to the kinetic 
energy of the translationally-localized linear-rotor impu¬ 
rity, with B the rotational constant and J the angu¬ 
lar momentum operator. In the absence of an external 
bath, the impurity eigenstates, |j, m), are labeled by the 
angular momentum, j, and its projection, m, onto the 
laboratory-frame z-axis. Unperturbed rotational states 
form (2j + l)-fold degenerate multiplets with energies 

Ej = Bj{j + i) [mHanoi- 



FIG. 1. Action of the canonical transformation, Eq. 
on the many-body system. Left: in the laboratory frame, 
(x,y,z), the molecular angular momentum, J, combines with 
the bath angular momentum, A, to form the total angular 
momentum of the system, L. Right: after the transformation, 
the bath degrees of freedom are transferred to the rotating 
frame of the molecule, As a result, the molecular 

angular momentum in the transformed space coincides with 
the total angular momentum of the system in the laboratory 
frame. 
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Many-body problems such as given by the Hamilto¬ 
nian 0 are typically hard to solve. The conventional 
approaches to tackle them include, when applicable, per¬ 
turbation theory, renormalization group, or in principle 
uncontrolled methods such as those based on the selective 
diagram resummations, as well as purely numerical tech¬ 
niques. An alternative, actively used since the develop¬ 
ment of classical mechanics, involves canonical transfor¬ 
mations of the underlying Hamiltonian [3Ql l3T] . Here the 
idea is to partially diagonalize the Hamiltonian and/or 
to expose the constants of motion, which allows to reveal 
some of the eigenstates’ properties exactly. In the con¬ 
text of impurity problems, typical approaches employ the 
collective bath variables as a generator of the symmetry 
transformations, as it has been used e.g. in the polaron 
theory [32ti3^ . 

In the angulon problem discussed in this paper, the to¬ 
tal angular momentum is a good quantum number. How¬ 
ever, due to the coupling of bath degrees of freedom with 
the impurity coordinates, as given by the third term of 
Eq. Q, this conservation law is not apparent. Here we 
introduce a canonical transformation which makes this 
constant of motion explicit and allows to achieve several 
other goals listed below. The corresponding operator S 
uses the composite angular momentum of the bath as a 
generator of rotation, which transfers the environment 
degrees of freedom into the frame co-rotating along with 
the quantum rotor. The transformation is given by: 


g — ^-iO^ky 


( 4 ) 


The angle operators, (0, ^, 7 ), act in the Hilbert space of 
the rotor, and 


^ = E (5) 

kXfiu 


well as the properties of the J' operator are presented in 
Appendix 

Let us now discuss the physical meaning of the trans¬ 
formation S. In order to describe the composite sys¬ 
tem, it is natural to introduce two coordinate frames, as 
schematically shown in Fig. The laboratory frame, 
(x,^, 2 :), is singled out by the collective state of the 
bosons, while the rotating impurity frame, {x'^y'^z')^ is 
defined by the instantaneous orientation of the molecular 
axes. The relative orientation of the two frames is given 
by the eigenvalues of the Euler angle operators, (0, ^, 7 ), 


acting in the impurity Hilbert space. The S operator 
transforms the many-body state of the bosons into the 
rotating molecular frame, using A as a generator of quan¬ 
tum rotations. In turn, as we show below, the molecular 
state in the transformed frame becomes an eigenstate of 
the total angular momentum of the system, which is a 
constant of motion. 

Introducing the body-fixed coordinate frame bound to 
the impurity makes explicit an additional quantum num¬ 
ber, n, which gives the projection of the angular momen¬ 
tum onto the rotor axis z'. The angular momentum basis 
states, |j, m,n), are therefore the eigen stat es of the J^, 
and J' operators, as given by Eqs. (C1)-(C3) of the 
Appendix. 

For a linear-rotor molecule in the absence of a bath 
the total angular momentum, L = J - 1 - A, coincides with 
J. Therefore, L is perpendicular to the molecular axis 
z', resulting in n = 0. With the bosons present, the 
total angular momentum is no longer perpendicular to 
z' ^ providing the molecular state with nonzero n in the 
transformed frame. In other words, the transformation 
0 converts a linear-rotor molecule into an effective ‘sym¬ 
metric top’ [ 20 ] by dressing it with a boson field. 

Compared to the original Hamiltonian, Eq. 0 , the 
transformed Hamiltonian, Eq. ^ , possesses the following 
properties: 


is the collective angular momentum operator of the 
many-body bath, acting in the Hilbert space of the 
bosons. Here cr^ denotes the vector of matrices fulfill¬ 
ing the angular momentum algebra in the representation 
of angular momentum A. 

The transformation brings the Hamiltonian 0 into 
the following form: 


H = S-^HS = B{3' - kf 


+ E ‘^ki>\\ffik\ti + y^y\{k) 

kXfj, kX 



+ bk\o 


(6) 


Here Vx{k) = Ux{k)x/{2X + l)/{4Tr) and J' is the 
‘anomalous’ angular momentum operator acting in the 
rotating frame of the impurity. Since the components of 
J' act in the body-fixed frame, they obey anomalous com¬ 
mutation relations uniES] as opposed to the ‘ordinary’ 
angular momentum operator, J of Eq. Q, which acts in 
the laboratory frame. The details of the derivation, as 


1. 'H is explicitly expressed through the total angular 
momentum, which is a constant of motion. Due to 
the isotropy of space, the eigenstates of the origi¬ 
nal Hamiltonian, H, are simultaneous eigenstates 
of the total angular momentum operators, and 
hz: and thus can be labeled as \L,M). The trans¬ 
formed states, S~^\L, M), are hence the eigenstates 
of the transformed Hamiltonian, H. As detailed 
in Appendix 0 these transformed states are also 
eigenstates of the J'^ operator with the eigenvalues 
L{L - 1 - 1), corresponding to the total angular mo¬ 
mentum. Consequently, the operator in Eq. 
can be replaced by the classical number 1 /( 1 / -h Ip 

2 . does not contain the impurity coordinates 
which allows to bypass the intractable angular mo¬ 
mentum algebra, arising from the impurity-bath 
coupling. The angle operators of the original 
Hamiltonian, Eq. Q, couple the impurity states 
with every single boson excitation, which results in 
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the problem of adding an infinite number of angular 
momenta in three dimensions. The latter involves 
working with Wigner 3nj-symbols of an arbitrarily 
large order. In the transformed Hamiltonian, on 
the other hand, the problem is reduced to adding 
the angular-momentum projections of the impurity 
and the bath. There, the impurity-bath coupling, 
J' • A, has the form of spin-orbit interaction and 
does not lead to an involved angular momentum 
algebra. 

3. i-L can he solved exactly in the limit of a slowly ro¬ 
tating impurity, H ^ 0, see Sec. m 

4. H allows to find the eigenstates eontaining an infi¬ 
nite number of phonon excitations, which is crucial 
e.g. to account for the maeroseopie deformation of 
the eondensate. This follows directly from sub. 
and is detailed in Sec. m 

5. H contains information about the deformation of 
the eondensate in the rotating impurity frame. 
Compared to the laboratory frame, where the de¬ 
formation of the bath is averaged over the angles, 
this provides an additional insight into the nature 
of the many-body state and, consequently, into 
the origin of the angulon instability, discussed in 
Sec, [ml 


III. MACROSCOPIC DEFORMATION OF THE 
BATH AND THE EMERGING INSTABILITY 



log[B/uo] 


FIG. 2. Change of the angulon spectral function, Al{w), 
where uj — E — BL{L + 1), with the rotational constant B, 
for three lowest total angular momentum states. The L > 0 
states show an instability in the spectrum. The red dashed 
line shows the deformation energy, Eq. (|^, which is inde¬ 
pendent of L. The circles indicate the points for which the 
phonon density modulation is shown in Fig. 


In the limit of a slowly-rotating impurity, B ^ 0, the 
Hamiltonian § can be solved exactly by means of an 
additional canonical transformation: 

= u-^nu (7) 


Therefore we introduce a variational ansatz based on 
single-phonon excitations on top of the bosonic state 
macroscopically deformed by the operator U : 

= gLMmiMQi)+ Y,akXnhixn\^)\LMn) ( 10 ) 

kXn 


where 


U = exp 


E 

. kX 


Vxjk) 

Wkx 



(8) 


with Wkx = cc;/c+HA(A+l). This transformation removes 
the terms linear in the bosonic operators, replacing them 
by the deformation energy of the bath. 


EA,i = -E^x{.kf/Wkx (9) 

kX 

As a consequence, in the limit of H = 0, the vacuum of 
phonon excitations, |0), becomes the exact ground state 
of Eq. 0. On the other hand, such a coherent shift trans¬ 
formation corresponds to a macroscopic deformation of 
the bath, and could not be easily performed on the orig¬ 
inal Hamiltonian 0 where the impurity coordinates are 
strongly coupled with the bath degrees of freedom. 

Here we are interested in the effect of a slowly-rotating 
impurity on the many-body state of the environment. 


The states of an isolated symmetric-top molecule are 
characterised by three quantum numbers: the angular 
momentum, L, its projection, M, onto the laboratory- 
frame z-axis, and its projection, n, onto the molecular 
symmetry axis, z'. For a linear rotor molecule, the an¬ 
gular momentum vector is always perpendicular to the 
molecular axis and therefore n is identically zero. The 
transformation Q, however, transfers the bosons to the 
molecular frame, thereby creating an effective ‘many- 
body symmetric-top’ state. The latter consists of a linear 
rotor impurity dressed by the field of bosons carrying fi¬ 
nite angular momentum. As a result, the total angular 
momentum of such a symmetric-top is no longer perpen¬ 
dicular to the linear rotor axis and provides the finite 
values of the projection n. See Appendix for more 
details. 

It is worth emphasizing that the non-transformed 
many-body wavefunction corresponding to Eq. (10) is 
given by \(l)) = S• U \fi). Therefore it is a highly-involved 
object with an infinite number of degrees of freedom en- 
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tangled with each other. The simple ansatz of Eq. (10) 
was made possible by the consecutive canonical transfor¬ 
mations, Eqs. 0 and ([^. Eurt her more, it is straight¬ 
forward to extend Eq. (10) to bath excitations of higher 
order, since this does not generate any complexities re¬ 
lated to the angular momentum algebra. 

Performing the variational solution for the energy, E = 
('0|^|'0)/('0|'0), we obtain the condition 


-EeBL{L^1)-El{E)=0 


( 11 ) 


which has the form of a Dyson equation with self-energy 
El{E) [36], as given by Eq. ( D13| ), see Appendix l^for a 
detailed derivation. Eq. (11) can be rewritten in terms of 
the angulon Green’s function as [Gl{E)]~ = 0, where 


[Gl{E)]-^ = [GI{E)]-^-El{E), (12) 

with [Gl{E)]-^ = -E + BL{L -f 1). 

The ground- and excited-state properties of the sys¬ 
tem are contained in the spectral function, Al{E) = 
lm[GL{E + '^50+)]. Without restricting the generality of 
what follows, we assume potentials whose angular mo¬ 
mentum expansion, Eq. §, is given by the Gaussian 
form-factors, fx{r) = and nonzero 

magnitudes, uq and ui, in two lowest angular momentum 
channels. We assume an anisotropy ratio of ui/uq = 5, 
a range vq = ri = 15 and set the interactions 

with A > 1 to zero. Eurt her more, we use a Bogo linbov- 
type dispersion relation, = x/^k{^k + where 

Ck = j2m with m the mass of a boson. We choose 
the boson-boson interaction ^bb = and 

density n = 0.014(mi4o)^/^. This choice of parame¬ 
ters reproduces the speed of sound in superfluid ^He 
for = 27r X 100 GHz |37|. Eig. shows the depen¬ 
dence of the spectral function on the rotational constant 
B for the three lowest rotational states. The width of 
the lines reflects the lifetimes of the corresponding lev¬ 
els. In Ref. [27|, we studied the non-transformed Hamil¬ 
tonian ^ using a variational ansatz based on single¬ 
phonon excitations. Using this ansatz we found that the 
angulon states become stable after crossing the phonon 
threshold at zero energy. Here this is no longer the case, 
since the transformation U of eq. <0 introduces an in¬ 
finite number of phonon excitations into the variational 
ansatz. This leads to a energetic renormalization of the 
phonon emission threshold providing all the excited an¬ 
gulon states with decay channels for phonon emission. 
This, in turn, leads to a finite lifetime for any magnitude 
of the impurity-bath coupling. 

In the limit of H ^ 0 the molecule is not rotating 
and is inducing an anisotropic deformation of the bath, 
corresponding to the mean-field-like deformation energy, 
Eq. The magnitude of the deformation energy de¬ 
creases with B monotonously and determines the gen¬ 
eral shape of the spectrum. Apart from the deforma¬ 
tion energy which is identical for all L’s, the energy of 
the angulon acquires an additional contribution due to 
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FIG. 3. Phonon density in the impurity frame for selected 
values of \og[B/ uq], which are specified in the right-top corner 
of the panels and (a)-(d) as labeled in Fig. (e) far to the 
right from the instability, at \og[B/uo] = 3.5 for L = 1 and 
at \og[B/uo] =3.0 for L — 2. The coordinates {x,z) are in 
units of ^ 
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phonon excitations in the surrounding medium. The lat¬ 
ter corresponds to the rotational Lamb shift discussed in 
Ref. |27], which has been observed as the renormalization 
of the rotational spectrum for molecules in superfluid he¬ 
lium nanodroplets [H]. Most importantly, we find that 
for the excited states with L > 0 there exists a critical ro¬ 
tational constant, where a discontinuity in the rotational 
spectrum occurs. This effect corresponds to a transfer of 
one quantum of angular momentum from the bath to the 
impurity. One can see that the faster the rotation (i.e. 
the larger L), the earlier this instability occurs. Such an 
instability has been briefly discussed in Ref. m, where 
it was referred to as Many-Body-Induced Fine Structure 
of the second kind. 

While the instability can be detected using spec¬ 
troscopy in the laboratory frame, an insight into its origin 
can be gained by making use of the canonical transforma¬ 
tion, Eq. Namely, in the frame co-rot at ing with the 
impurity, the instability manifests itself as a change of the 
phonon density, (6j6r); for analytic expressions see Ap¬ 
pendix Fig. shows the phonon density for L = 1 and 
2 at five different values of the impurity rotational con¬ 
stant. Darker shade corresponds to higher density. Far to 
the left of the instability, panels (a), the impurity is rotat¬ 
ing slowly and the bosons are able to adiabatically follow 
its motion. As a result, the surrounding bath becomes 
polarized, which manifests itself in a highly-asymmetric 
phonon density. The shape of the density modulation is 
given by the first spherical harmonic which arises due to 
the A = 1 term in the impurity-boson potential included 
in our model. Closer to the instability, panels (b), the 
phonon density increases, signaling the onset of the reso¬ 
nant phonon excitations. At the right edge of the insta¬ 
bility, panels (c), the phonon density drops drastically. 
Further away from the instability, the density distribu¬ 
tion becomes the more symmetric the faster the impurity 
rotates, as illustrated in panels (d), (e). In other words, 
when the rotational constant exceeds the critical value 
given by the instability, it becomes energetically unfavor¬ 
able for the bosons to follow the motion of the impurity. 
As a consequence, the bosonic bath does not possess fi¬ 
nite angular momentum, which results in the spherically- 
symmetric density distribution. Thus, the phonon den¬ 
sity in the transformed frame can serve as a fingerprint 
of the angular momentum transfer from the bath to the 
impurity which takes place at the instability point. 

It is important to note that the ‘angulon instability’ 
discussed here is fundamentally different from the vor¬ 
tex instability [29], also associated with rotation. The 
comparison between the two is summarized in Table [T| 
First, the rotation of the impurity is inherently three- 
dimensional and does not involve any specific rotation 
axis. This is different for a vortex line which singles out 
a particular direction in space. Second, the formation of 
a vortex requires a transfer of one unit of angular momen¬ 
tum per particle in the bath. In the angulon instability, 
on the other hand, a finite (small) number of rotational 
quanta is shared between the impurity and the collec¬ 


TABLE I. Comparison of the angulon instability with the 
vortex instability. 



Angulon 

Vortex 

Corresponding rotation 
Angular momentum transfer 
Circulation 

Spherical, L^ 
h 

zero 

Planar, L^ 
h per particle 
integer x27rd/?Ti 


tive state of the many-particle environment. Finally, the 
vortex instability leads to a finite circulation around the 
vortex line, which is absent for the angulon instability. 


IV. EXPERIMENTAL IMPLEMENTATION 


The described effects can be observed experimentally 
both with molecules trapped in strongly-interacting su¬ 
perfluids, such as helium droplets m, or molecular im¬ 
purities immersed in weakly-interact ing Bose-Einstein 
condensates [29]. The dependence of the angulon self¬ 
energy, of Eq. (12), on the many-body parameters 


can be revealed by measuring the relative shift between 
the rotational states of a diatomic molecule. Since the 
effects will be most pronounced for the molecular states 
possessing a small rotational constant 5, experiments 
involving molecules in highly-excited vibrational states 
provide the most natural setup. In the context of ul¬ 
tracold gases, the latter include photoassociation spec¬ 
troscopy [38] and measuring nonzero angular momentum 
Feshbach resonances m- In both cases, the shifts of 
the spectroscopic lines will be proportional to the an¬ 
gulon self-energy, as schematically illustrated in Fig. 

An alternative possibility is measuring as a shift 
of the microwave lines in the spectra of weakly bound 
molecules [40], prepared using one of these techniques. 
In frequency domain, at sufficiently low temperatures the 
width of the lines will correspond to the angulon lifetime. 
The instability shown in Fig. corresponds to the van¬ 
ishing quasiparticle weight with a related emergence of 
a broad incoherent background and therefore can be de¬ 
tected as a line broadening with increasing impurity-bath 
interactions. In the time domain, on the other hand, the 
angulon Green’s function can be detected using Ramsey 




FIG. 4. Detection of the angulon self-energy Sl using 
(a) photoassociation spectroscopy [38] and (b) shift of p- and 
d-wave Feshbach resonances ES]. 
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and spin-echo techniques mnn]. In such measurement, 
the angulon instability will leads to dephasing dynam¬ 
ics with a related pronounced decay of the Ramsey and 
spin-echo contrast [msi]. 

While in superfluid helium the interactions cannot be 
tuned as easily as in ultracold gases, the range of chemical 
species amenable to trapping is essentially unlimited | 21 | . 
The latter, combined with the advances in the theory of 
molecule-helium interactions [43] paves the way to study¬ 
ing angulon physics in a broad range of parameters. 


V. CONCLUSIONS 


In this paper, we studied the redistribution of orbital 
angular momentum between a quantum impurity and a 
many-particle environment. We introduced a technique 
which allows to drastically simplify the problem of adding 
an infinite number of angular momenta which occurs in 
the regime of strong interactions. The essence of the 
method - a novel canonical transformation - paves the 
way to eliminating the complex angular momentum al¬ 
gebra from the problem, as well as to exposing the prob¬ 
lem’s constants of motion. We exemplified the tech¬ 
nique’s capacity by studying an instability which occurs 
in the spectrum of the many-particle system due the in¬ 
teraction between the bath and the rotating impurity. 
Such an instability should be detectable with molecules 
in superfluid helium droplets m and might be responsi¬ 
ble for the long timescales emerging in molecular rotation 
dynamics in the presence of an environment [44], which 
presently lacks even a qualitative explanation. More¬ 
over, the rotating impurities can be prepared experimen¬ 
tally in perfectly controllable settings, based on ultracold 
molecules immersed into a Bose or Fermi gas [laiiHiss] 
and cold molecular ions inside Coulomb crystals [23]. It is 
important to note that the transformation, as defined by 
Eq. Q, is quite general, and can be applied to extended 
Frohlich Hamiltonians [ 36 ], to impurities with complex 
rotational structure [20], Rydberg molecules [^1461148] . 
as well as to the case of a Fermionic bath m- 

The ultimate goal of our approach is to find a series of 
canonical transformations that would lead to exact solu¬ 
tions to the many-body Hamiltonians of the same class as 
Eq. §. This resonates with Wegner’s idea of the contin¬ 
uous unitary transformations m, which underlies one 
of the Hamiltonian formulations of the renormalization 
group approach [51]. 

Finally, the impurity problem considered here can be 
used as a building block of a general theory describing the 
redistribution of orbital angular momentum in quantum 
many-particle systems. This opens up a perspective of 
applying the techniques of this article to the several prob¬ 
lems in condensed matter m and chemical m physics. 
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Appendix A: The angular momentum representation 


The creation and annihilation operators of Eq. 0 
are expressed in the angular momentum representation, 
which is related to the Cartesian representation as: 

(Al) 

4 = (A2) 

Xfj, 


The quantum numbers A and fi define, respectively, the 
angular momentum of the bosonic excitation and its 
proj ectio n onto the laboratory-frame z-axis. Eqs. 
and (A 2 ) correspond to the following commutation rela¬ 


tions: 


[6k,&L] = (27r)3<5(3)(k-k') (A3) 


[^fcA/ii 4'A';<'] “)<^aaA^^' (A4) 

In the coordinate space, the transformation between 
the representations is defined as: 

^ ^r) (A5) 


(A6) 

A/x 

with the corresponding commutation relations: 

[&r, - r') (A7) 

~ ^ )^XX'Sf^fi' (AS) 

The operators in the coordinate and momentum space 
are related through the Eourier transform. 



from which one can obtain the corresponding relation for 
the angular momentum components 


bl=i\- 




j kdkjxikr) 


with j\{kr) the spherical Bessel function |S3]. 


(AlO) 
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Appendix B: The canonical transformation 


where i^k = { — 1 , 0 ,+!}, and with the rotation opera¬ 
tors: 


Here we provide details on the derivation of the trans¬ 
formed Hamiltonian, Eq. §). 

In the angular momentum representation, the boson 
creation and annihilation operators, bl.^^ and are 

defined as irreducible tensors of rank A [28]. Conse¬ 

quently, they are transformed by the ^-operator of Eq. 
@ in the following way: 

S-%x,S = (Bl) 

V 

= (B2) 


= (-1)'=+! + e, 7 ) (B 8 ) 


Jk,D^M8n) 


Here Ci 




= v'A(A + i)c^;;;y,i?^;,,,(<^,0,7) (b9) 

^ .1 ^ are the Clebsch-Gordan coefficients 1281. 
By usin g th e latt er property, one can show that the 
operators (B4)-(B 6 ) transform under Eq. Q in the fol¬ 
lowing way: 


Here ^, 7 ) are Wigner D-matrices [28| whose ar¬ 

guments are the angle operators defining the relative ori¬ 
entation of the impurity frame with respect to the labora¬ 
tory frame. These expressions can also be derived using 
the explicit expression for the angular momentum of the 
bosons, Eq. (§. 

The Wigner rotation matrix appearing in Eq. pri is 
complex conjugate with respect to the one of Eq. (B2) 
and therefore corresponds to an inverse rotation. As a 
result. 


(E = E (B3) 

/X fl 

and the second term of Eq. § does not change under 
the transformation. 

Similarly, in the last term of Eq. 0 we use that 
Ea/x (^5 0 ) = 0^ 0 ), which leads to cancel¬ 

lation of the Wigner E)-matrices. In such a way, the 
transformation S eliminates the molecular angle variables 
from the Hamiltonian. 

The transformation of the molecular rotational Hamil¬ 
tonian, HJ^, turns out to be slightly more cumbersome. 
In the laboratory frame, the angular momentum vector is 
defined by its spherical components, J = {J_i, '>^+ 1)5 
where: 


Ji = S-^JiS = Ji- (BIO) 

/c=-l, 0 ,l 

After some angular momentum algebra, we obtain the 
following expression for the square of the angular mo¬ 
mentum in the transformed frame: 

S-^3^S = Ji - J+i - J-iJ+i = (J' - kf (Bll) 

Here J' is the angular momentum operator in the rotat¬ 
ing molecular (i.e. body-fixed) coordinate frame [35] . 
which can be expressed via the laboratory-frame compo¬ 
nents as: 

(B 12 ) 

k 


The spherical components of J' are expressed through 
the Cart esian components using the relations analogous 
to Eqs. (B4)-(B6). Note that this makes the J' oper¬ 


ators different from the so-called contravariant angular 
momentum components used by Varshalovich |28| . 

The molecular-frame angular momentum operators 
obey the anomalous commutation relations with one an¬ 
other [T^ [ 35 ] . 


J' f' 




(B13) 


and the following commutation relations with the rota¬ 
tion matrices: 


Jo = J 2 

j+i = - 


j-i = 


V 2 ^ J) 

Jx ^Jy^ 


^/2 ( 


(B4) 

(B5) 

(B6) 




= - v/a(a + i) 7 ; 7 ,d ^,,+,{$, e, 7 ) {B14) 




see Refs. [23 [35]. We use the analogous notation for 
the components of the total angular momentum of the 
boso n s A = {A_i, Aq, A+i}, Eq. The operators 

( |B4[ )-( [B^ obey the following commutation relations with 
each other: 


Jii Jk 




(B7) 


= (- 1 ) VA(A + 0, 7 ) {B15) 

It is worth noting that in the case of a linear-rotor 
molecule, the molecule-boson interaction does not de¬ 
pend on the third Euler angle, 7 . However, this angle 
must be preserved in Eq. Q , as well as in all the deriva¬ 
tions described above, in order to keep the transforma¬ 
tion unitary. 






























Appendix C: Molecular states in the transformed 
space 

In the main text and Fig. we have introduced two 
coordinate frames: the laboratory one, {x^y^z)^ and the 
molecular one, A general molecular state, 

therefore, can be characterised by three quantum num¬ 
bers: the magnitude of angular momentum, j; its pro¬ 
jection, m, onto the laboratory-frame 2 :-axis; and its pro¬ 
jection, n, onto the molecular-frame z'-axis: 

J2|j,m,n) =j(i +1)1 (Cl) 

Jz \j, m, n) = m\j, m, n) (C2) 

J'z \j, m, n) = n|j, m, n) (C3) 

In the angular representation, the corresponding wave 
functions are given by m 


{4>,e,-i\j,m,n) = i4>, 0, 7 ) (C4) 

The action of the space-fixed and molecule-fixed com¬ 
ponents of angular momentum is given by the general 
formula nms]: 

Jk \j, m, n) = vWTl)C'j;™.+^ I j, m + k,n) (C5) 

Jk\h m, n) = (-I)''+ l)C'j;^;T,-fcb) rn,n-k) 

(C 6 ) 

where k = { — 1,0,+!}. Thus, in the molecular frame the 
raising operators lower the projection quantum number 
n and the lowering operators raise it. 

Unlike for nonlinear polyatomic molecules [20], the an¬ 
gular momentum of a linear rotor is always perpendicular 
to the internuclear axis (defining z')^ and therefore n is 
identically zero. However, this is the case only before the 
transformation S is applied. Let us consider the most 
general many-body state in the non-transformed frame, 

\L, M) = J2 4xjC^:r^,xJjm0) ® \kXfi)i (C7) 

kXfi 

jm]i 


The molecular states \jmQ) are the eigenstates of the 
mole c ular angular momentum operator, as given by Eqs. 
(C1)-(C2). The same relations are fulfilled for the col¬ 
lective bosonic states: A^|/cA/i) = A(A + l)\kXfi) and 
Az\kXfi) = y\kXfi), where A is defined by Eq. (^, and k 
is the linear momentum. The index i labels all the pos¬ 
sible boson configurations resulting in a collective state 
|/cA/i), spanning the complete many-body Hilbert space 

of the bosonic bath. _ 

It is straightforward to show that the state (C7) is 
an eigenstate of the total angular momentum operator, 
L = J + A: 


(C8) 

(C9) 


By acting on jL, M) with after some angular mo¬ 
mentum algebra, we obtain the state in the transformed 
frame: 

S-^\L,M) = fixnlLMn) ® \kXn)i (CIO) 

kXni 

where the coefficients are given by 

fkXn = ^ ^kXj^L^-n]X,n (Cll) 

j 

We see that the transformation effectively transferred 
the angular momentum of the bosons to the molecular 
frame. This is reflected by the fact that the transformed 
state, becomes an eigenstate of the body- 

fixed angular momentum operator, J'^, with the eigen¬ 
values of the total angular momentum operator, L^, i.e. 

J'^ =L{L + 1) (C12) 

Each state \LMn) in the superposition of Eq. ( |C10[ ) is 
an effective symmetric-top state [20], with the projection 
of total angular momentum on the molecular axis entirely 
determined by the boson field. 


Appendix D: Derivation of the Dyson equation from 
the variational principle 


We minimize the energy obtained from the expectation 
value of Eq. 0 with respect to the variational state: 

li’) = gLMmiMO) + Y,(^kXnbiMLMn) (Dl) 

kXn 

Minimization with respect to OL%xn 9lm yields the 
following equations: 


[—E -f BL{L -1- 1)] qlm + BxJL{L + ^kXO^kXn = 0 

kX 

(D2) 

and 


[—E BL{L + 1) + Wkx] GLkXn — ^ nv'n nv^kXu 

V 


EBSn,±l^kX y^^^k'X'O^k'X'n 

k'X' 

= —BxJ L{L + l)^kXgLM^n,±l 


(D3) 


where we defined 6n,±i 
x/X{X^l)Vx{k)/Wkx. ^Wkx = ojk 


— ^n,l T ^n, —15 C/cA — 
■h HA(A "T 1), and 


-.L _ 


= {L Mn\ J ' \LM u). In what follows, we show that 


Eqs. ( |D2[ ) and ( |D3[ ) can be solved in closed form. 

Ei rst, the angular-momentum coupling term of 
Eq. (D3) is given by: 


+ i\/A(A + 1) - ix(ix + l) 

X y L{L + 1) — + l)(5„^y_|_i 

+ ^^yx{x + ^)-u{u-l) 

X ^LiL + l)-uiu-l)Sn,.-i (D4) 


V\L,M) = L{L + 1)\L,M) 
Lz\L,M) = M\L,M) 
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Assuming that Vx{k) 7^ 0 for A = 0,1 only, we obtain 
that Eqs. (D2) and (D3) are solved by akXn = 0 for 


A = 0. Consequently, and akio are the only 

variational parameters. 


For akio we obtain 

[—E + BL(^L + 1) + LUk + 2B] cxkio 


—Bx/2L{L + l)((a/cii + OLki-i) — 0 


(D5) 


For the aki±i components we find two identical equa¬ 
tions 

[—E + BL{L -f 1) + cj/c] 

—Bx/2L{L -1- l)(a/cio + B^j^i E 


= - 5 a / 1 /( 1 / + i)^kigLM 


(D6) 


By symmetry we expect \akii\ = |(a/ci-i|, however, if 
akii = —were true, Eq. ( |D5| ) would imply a/cm = 
0. This in turn would lead to a contradiction in Eq. (D6) 
which shows that akii = cxki-i- 
Thus, from Eq. (|D5|) we obtain 


<a/cio = 


2B^2L{L +1) 




—E T uj]^ T BL(^L -f-1) T 2B 
Let us now define the inverse propagator 

AB‘^L{L +1) 


(D7) 


Psik) — i?-L(-Zj+l)—-B+Wfe- 


—E -\- LOj^ BL(^L “1“ 1) “1“ 2B 

(D8) 


and rewrite Eq. (D6) as: 

BCk 


CTfell — — 


PE{k) 


1 * P' 1 /L{L + 

^ 2_^ - Tr7Tr\ - 9lm 


PE{k) 


(D9) 

In addition it is convenient to introduce the variable y as 
9lmX = X] ^kioikii (DIO) 


After multiplying Eq. (D9) with and integration over 
k we find 


= -B'/L{L + l) 


j^dke^jPEik) 

l + B!^dkilJPE{k) 


(Dll) 


This finally yields the Dyson equation 

- 5 + BL{L + 1) - El(5) = 0 (D12) 

where the self-energy is given by 

///dkilJPEik) 


Ei(E) =B"L(L+1) 


l + Bj^dkekJPE{k) 


(D13) 


and 


^ki = V2 


vm 

cj/c + 2B 


(D14) 


Furthermore, we absorbed the deformation energy Edef, 
Eq. ([^, which is identical for all the L-levels, into the 
definition of E. Note, that if B/ux ^ 1, the self energy 
El -> BL{L 1) and the Dyson equation is solved 
by 5 = 0. This means that for weak interactions the 
impurity levels are shifted by the mean-field deformation 
energy only. 


The self-energy of Eq. ( |D13[ ) can be partially evaluated 

(D15) 


analytically. It is convenient to define 
uj = E - BL{L + 1) 


and to rewrite the retarded self-energy, E5f^(u;) = El{oj-\- 
iO^), as 


= 2B^L{L 


.N Xl{u}) 

'^1 + 2xl(w) 


where 


Xl(w) 


-f 


dk 


Vi{kf 


[ujk + 25]2 ^cD+io+(^) 


(D16) 


(D17) 


The integrand of Xl{^) possesses poles at the momenta 
ko satisfying ujkQ = u; for L = 0 and at the momenta ki ^2 
satisfying ujj^^ = cj + 2BL and ujk 2 = oj — 2BL{L + 1) 
for states with L > 0. Using the relation l/{x = 

ViXjx) — i7rS(x) this reveals the onset of the scattering 
continua in the spectral function. 

For L = 0 one finds 


where 


Co = 


ImXL=o(w) = 7r6>(w)Co 

Vi(ko/ 


[w + 2B]2 

while for L > 0 one has 

TT 

ImXL>o(w) = -d(aJki) 




'dtok 

-1 

_ dk _ 



k=ko 


1 


y/l + 4D(iy + 1) 
1 

yi + 4L(L +ly 


where 


Cl,2 — 


U(fco)2 

■2Bf 


\U} ■ 


~dL0k 

-1 

. dk _ 



(D18) 

(D19) 

Cl 

C2 

(D20) 

(D21) 


k=k\ 


Finally, the real part of Xl{^) follows from the principal 
value integration. 
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Appendix E: Deformation of the phonon density 


From Eq. (A6) we obtain the expression for the phonon 
density in the rotating impurity frame: 


n(r) = (blbr) 

Xfi 

X' fj,' 

(El) 


Using Eq. 
butions: 


(AlO), we evaluate the partial-wave contri- 


= j kdk J k'dk' jx{kr)jx'{k'r) {blx^h'X'n') 

(E2) 


We calculate the expectation values, (...), with re¬ 
spect to the states in the transformed frame, \(Q_ = U\ 2 p)^ 
where U and \ 2 p) are given by Eqs. ^ and (10) of the 
main text. Einally, the expectation values of (b\xJ^k'X'fi') 
are given by: 


iH:XfJ^k' X' ^d') — 


Vx{k)Vy{k') 

Wkx Wk^x' 


~9LM^k'X'0 


Vx{k) 

Wkx 


- QlmoiIxo 


Vx'{k')Vx{k) 
Wk'X' Wkx 


+ WkXii\^ 


(E3) 
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